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Abstract 

A chemical system consisting of two species, one of which evolves deterministically and indepen- 
dently of the other, which in turn is driven by the dynamics of the former and by an additional 
multiplicative Gaussian white noise, displays a 1// noise for intermediate to large frequencies. 
A novel mechanism responsible for the 1/f noise is suggested. 

PACS numbers: 05.40. Ca 
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Although systems exhibiting 1//" noise, and related systems which display power-law 
type distributions, are abundant in nature, there is no generally accepted explanation for 
their origin and ubiquity. Self-organized criticality has been proposed as such an expla- 
nation, but it has been pointed out that the notion of self-organized criticality and 1//" are 
mutually exclusive in many cases 0. Recently the emergence of power-law distributions and 
long-time temporal correlationshas have been explained for example in terms of a general- 
ized logistic systems P]. More importantly, Davidsen and Schuster in Ref. @] have proposed 
a mechanism explaining the origin of 1//°, a ~ 1, noise in the low-frequency regime, which 
corresponds to very long lasting temporal correlations. The systems considered in Ref. Q 
were driven by a Gaussian white noise (GWN) — in other words, no preexisting power-law 
distributions or fractal noises have been assumed, even implicitly. Power-law distributions 
have also been shown to result, at the level of the mathematical formalism used, from non- 
linear transformation between various stochastic variables for a very wide class of underlying 
"fundamental" distributions 0. 

In the present Letter we consider a simple chemical system consisting of two species. Dy- 
namics of one of the species. A, is subjected to a multiplicative GWN and is further affected 
by the dynamics of the other species, B. The dynamics of B is purely deterministic and is 
not influenced back by the dynamics of A. Such a system can be realized experimentally. 
We show that our model system exhibits 1// noise for intemediate to large frequencies. 

Consider two chemical species, A and B. Let B decay autocatalytically and catalyse a 
decay of A: 



A + B — ^ A*i +B, (1) 



K 



B + B ^ B*i +B. (2) 

kb 

The downarrows mean that species A*, B* are stable and do not enter any reactions of 
interest. If the amounts of reagents are periodically incremented from outside 



Ja{t) = aoJ2 - , = feo 5^ 5(t - nT) , (3) 

n n 

the kinetics of reactions (|I])-(0) can be expressed as 
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a = -Kah + ]a{t) , h = -hb^ + jb{t) , (4) 

where a, b are concentrations of A and B, respectively. We further assume that the rate at 
which A decays fluctuates around a deterministic value 



K = ka + nriit) , (5) 

where rjit) is a GWN with {ri(t)) = 0, {rj{t)ri{t')) = a'^6(t — t'), and all higher correlations 
factorize. Reactions (|D-(@) with the above constrains can be realized experimentally in a 
flow reactor, for instance by a photoactivated chemical reaction with a radiation field acting 
as a source of noise 0; see also [0 and references quoted therein. 

Because of the 5-terms in the fluxes, we solve in a stroboscopic representation. Be- 
tween the pulses 



+ nT<t<(„ + l)T, (6) 

where 6„ = b{t = nT^). The pulses simply increment the concentrations by ^^nd 
respectively. Thus 



t 



nT+ 

nT <t<{n + l)T, (7b) 

where a„ = a(t = nT~^), fi = ka/ki, is a ratio of deterministic reaction rates and 

(n+i)r- 

E„ = exp I - y Kri{t')b{t')dt' I . (8) 



The {bn} form a deterministic series convergent to 



Note that > bo- Contrariwise, the series {a„} is stochastic and its reahzations need not to 
be convergent; only the series of expectation values {(a„)} can. Both a„ and En are random 
numbers, but since a„ depends only on times prior to nT, an and En are defined on disjoint 
intervals and for 77 (t) beging a GWN, the expectation value of their product factorizes. We 
have 



{En) = exp ( |kV^ 



1 + bnkfjT 



(10) 



(cf. 



and 



In the limit n — > 00 we replace 6„ by 600 and ([TTl) yields 

(12) 

There is a maximal addmissible noise level 



1- (1 + 600 exp 



n'ai^^=-^Hl + b^k,T) (13) 

beyond which the process diverges. 

We also need to specify a variance of the process a„: 



(El) + {anf {{ED - (Enf) 



(14) 



(1 + 6„A;feT)2M 

It is easy to verify that the series {f^} converges for k^ct^ < \i^^cr'^i^_^ and diverges for higher 
noise levels. 
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FIG. 1: Power spectra corresponding to the correlation function ([l9|): (i) ;U = 1, 60 = 1, (ii) = 2, 
60 = 1, (iii) A* = 2, 60 = 10, (iv) = 5, 5o = 10. Other parameters are ao = 1, = 1, T = 4. In 
each case k^o"^ = 0A5 n'^a'^^^. Note a changing slope of the spectrum (iv). 

Finally we need to examine the autocorrelation structure of the process a{t). We define 



C{t, r) = {a{t)a{t + r)) - (a(t)) {a{t + r)) . (15) 

For a GWN we obtain 

m— 1 

C{nT, niT) = {{al) - (a„)') J] (E„+,) , (16) 

i=0 

where Bi = + Note that all Ei^s in (|16D are defined on mutually disjoint intervals 

and depend on times later than nT. For n ^ 1, when the process reaches its stationary 
state, ( [I6| ) simplifies to 

provided that k^ct^ < ^k^ct^j^^. 

The Fourier transform of the autocorrelation function ( p!7D is by the Wiener-Khinchin 
theorem related to the power spectrum of the process {a„}, n ^ 1, averaged over realizations 
of the noise: 
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FIG. 2: Time-dependent spectrum of the system not sustained from outside. For all times t, the 
power spectrum displays a clear power-law behavior for large oj. Parameters are /U = 1, 60 = Ij 
ao = 1, feb = 1, K^cr^ = 1.98. Units oi uj, t are reciprocal, but otherwise arbitrary. 



E 

m=0 



1 — a cos uoT 



(18) 



1 — 2a cos uoT + ' 

where a = (1 + bookbT)~f^ exp {h^K^a'^ j {2kb)). For all noise levels such that is finite the 
series in ( ]T8| ) is convergent. Because ( pTSD corresponds to the power spectrum of a discrete 
time series "sampled" with a time step of T, there is a finite Nyquist frequency tt/T. If 
we want to probe higher frequencies, we need to examine correlations in the process a{t) 
between the infiuxes of the reagents. It is sufficient to calculate C{nT, r) with r = rriT + r*, 
where m is natural and < r* < T. For n ^ 1, C{nT, r) ~ C{t) with 



C(r) = ^ exp ( Wa"" \ . (19) 

The correlation function ( p!9|) corresponds to the power spectrum of a{t) averaged over 
realizations of the noise. Typical spectra are plotted in Fig. |l[ For low frequencies, uT < vr, 
which correspond to large r, the spectra are dominated by the behavior similar to (0). For 
very large frequencies (very small r), all curves plotted display a clear 1// (or l/oj in our 
notation) decay. This effect is generic: If we expand ([T9|) in powers of , we get 
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C{t) ^ Crn [1 - hU^^h - + O ((r^)) . (20) 

For very small r*, the correlation function approaches a non-zero value, which is responsible 
for the 1// decay at very large frequencies. The second term in the expansion would lead 
to 1/ P decay for intermediate frequencies; however, as r* cannot grow too large, this effect 
can show up only if the coefficient is large enough (cf. spectrum (iv) on Fig. |l]). If the time 
interval, T, between the consecutive influxes is large, r* can also be large and (|1^) simplifies 
to 



The coefficient in the above expression can be made sufficiently large by choosing appropriate 
values of the parameters of the system. 

The spectra of Fig. become flat for very high frequencies, but this is a numerical effect 
resulting from the roundoff errors acting like a white noise. 

Let us compare the above results with the situation in which the reactions are 
not sustained from outside, ja = jb = ^- Such reactions can be experimentally realized in a 
closed reactor, as opposed to the flow reactor considered above, but formally this corresponds 
to taking T ^ 1 and considering the flrst interval only. Under these assumptions we obtain 



{a{t)) = Z{t)e'^'\ (22a) 

(22b) 

where 



2 'l + hQhx){l + hQh{x+y))' 



If r = 0, the second of equations (p2D gives the variance of the process a{t). We can see 
that for small t the variance grows linearly, then reaches a maximum, and then decays; 



for long times this decay has a power tail ~ regardless of the noise level: for times 
large enough almost all realizations of the process a{t) vanish, and therefore a statistical 
difference between these realizations vanishes as well. This is clearly related to the presence 
and behavior of b{t). Recall that a decrease in b{t) limits the rate at which the substance A 
decays. For long times, when b(t) ~ 0, even wide fluctuations of the reaction rate K have 
little effect on the actual decay of A. If the concentration of B were kept constant {kb = 0), 
the expectation value of a{t), its variance and correlations described by (l2^) would either 
exponentially go to zero or expenentially diverge, depending on the noise level. 

Because in the not sustained case the process a{t) does not reach any stationary state, the 
correlation function ( P2b|) depends on two arguments and the corresponding power spectrum, 
P{uj,t), is also time-dependent. In ( P2b| ) we easily recognize the same type of small r 
behavior as in (^) which leads to the 1// tail for large frequencies. Indeed, this is what we 
observe numerically (Fig. |^). Note that this power tail does not change after averaging over 
the time, t: 



T 

1 



PH = - J P{uj,t)dt (23) 



with T ^ but finite, also displays this type of large frequencies behavior. We may also 
expect a slope of 1//^ type for intermediate frequencies and the parameters of the model 
large enough, but as time, t, increases, the frequency range corresponding to such a slope 
decreases, and for large (actually, not even very large) values of t the effect of the second 
order in r becomes negligible (Fig. For very large values of the parameters and small 
t, even the l/p decay appears for intermediate frequencies, but this effect also disappears 
as time, t, grows. The presence of the 1//^+" slope is, therefore, at most a transient effect, 
which does not survive time averaging (^Bf). In the limit t oo, the correlation function 
(p2b|) becomes 



which bears a formal similarity to (^TJ), but as the coefficient is very small, the r-dependence 
becomes practically unobservable. As we have said, this results from the fact that almost all 
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FIG. 3: Time-dependent spectra of the system not sustained from outside. Parameters are: (a) fx = 
l^bo = l,t = 1/2048, (b) /i = 5, 6o = 10, t = 1/2048, (c) ^ = 5, 6o = 10, t = 1/256, (d) ^ = 5, 
bf) = 50, t = 1/1024. Other parameters, common for all the presented spectra, are: ao = 1, kb = 1, 
n'^a'^ = 1.98 (except for (d), where K^a'^ = 0.5). Parameters of the spectra (b), (c) are the same as 
those of the spectrum (iv) of Fig. |^. Units of t are reciprocal, but otherwise arbitrary. 

realizations of the process a{t) vanish in the limit t ^ oo if the reactions are not sustained 
from outside. 

If the concentration of B is kept constant and the concentration of A is not incremented 
from outside, 

C{t, r) = ao exp {-2b{ka-bK^a^)t - b{ka~bK^a^ /2)t) . (25) 

A similar expression can be found if the concentration of A is incremented from outside. 
There is no power-law dependence on r in (^). It is now clear that the power-laws in ( pT]) 
and ( |2^ ) result from the fact that the stochastic process a{t) is driven by an "external", 
deterministic process b(t). 

In this Letter we have considered two cases: when the reactions are sustained from outside 
by periodic influxes of the reagents and when they are not. In the former case, the effects 
of noise are partially masked by the influxes, but since the noise effetcs in the latter cannot 
be observed due to the eventual decay of the reagents (Eqns. (0)), only the former admits 
a practical experimental realization. The mechanism responsible for the emergence of the 
1// noise in the system under consideration is following: kinetics of a proces driven by a 



multiplicative GWN (the decay of A) is additionally driven by another, deterministic process 
(the decay of B). To our knowledge, this is the first model in which a clear physical cause of 
the 1// noise has been established. As no power scalings or other special features have been 
assumed, even implicitly, our results suggest that this mechanism can have a more universal 
character. This idea will be pursued in a future research. 
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